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Abstract. The relation between quasar variability and 
parameters such as luminosity and redshift has been a 
matter of hot debate over the last few years with many 
papers on the subject. Any correlations which can be es- 
tablished will have a profound effect on models of quasar 
structure and evolution. The sample of quasars with red- 
shifts in ESO/SERC field 287 contains over 600 quasars 
in the range < z < 3.5 and is now large enough to 
bin in luminosity and redshift, and give definitive mea- 
sures of the correlations. We find no significant correla- 
tion between amplitude and redshift, except perhaps at 
very low redshift, but an inverse correlation between am- 
plitude and luminosity. This is examined in the context of 
various models for quasar variability. 
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1. Introduction 



One of the most important constraints on the structure 
of quasars is variability. Short term variations set limits 
on the size of the emitting region, and differences in 
the nature of the variation in the Xray, optical and 
radio domains give clues to the underlying structure. 
Variations on longer timescales of a few years are not so 
easily accounted for by the current black hole paradigm, 
and microlensing has been put forward as an alternative 
explanation for the observations. Although the most 
extensive monitoring of quasar variability has been done 
in optical wavebands, there is still little consensus even 
over the broad picture. One of the first problems is to 
parametrise the variation in such a way that it can be 
compared with models, or even more general expecta- 
tions. The timescale of variation and the amplitude are 
the two parameters which have been mostly studied, 
although some authors have succeeded in confusing the 
two. This typically involves claiming that in a short 
run of data, objects varying on a short timescale will 
achieve a larger amplitude than those varying on a longer 



timescale, and so amplitude can be taken as an (inverse) 
measure of timescale (Hook et al. 1994). 



In this paper we shall confine our attention to the 
correlation of amplitude with other parameters. Am- 
plitude is an easy parameter to measure, and involves 
none of the difficulties inherent in estimating timescale of 
variation. There are modes of variability in which some 
uncertainty will be introduced by the length of the run of 
observations, but this is a problem which can arise in any 
time series analysis. Most attention has so far been given 
to possible correlations of ampli tude with redsh ift or 
luminosity ( eg Tr cvese et al 
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Giallongo et al. 



Cristiani et al. 1996 ), but even in such 
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straightforward situations there has been little agreement. 

The reasons for the lack of consensus are not easy 
to understand, but it is clear that any effect which is 
present is not large compared with the cosmic scatter in 
the data, which in the case of amplitude appears to be 
about one magnitude, much larger than any photometric 
errors. Perhaps the most likely cause of disagreement is 
in the selection of samples for analysis. In order to obtain 
a meaningful measure of correlation it is essential that 
unbiassed samples are used, and that they cover a large 
range of redshift and luminosity. In this paper particular 
attention is paid to the sample selection process, and by 
extending the analysis to high redshifts the baseline for 
the measurement of correlations is greatly extended over 
earlier samples. 



2. Quasar samples 

2.1. Observational material 

The parent sample of quasars for the analysis of variabil- 
ity amplitude comes from a large scale survey and mon- 
itoring programme being undertaken in the ESO/SERC 
field 287, at 21h 28m, -45°. Quasars have been selected 
according to a number of criteria, including colour, vari- 
ability, radio emission and objective prism spectra, or a 
combination of these techniques. Redshifts for over 600 
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quasars have so far been obtained, and several complete 
samples defined within specific limits of magnitude, red- 
shift and position on the sky (Hawkins & Veron 1995| ). 
A detailed description of the survey is given by Hawkins 
& Veron (1995) and Hawkins ( 1996| ). Briefly, a large set 
of UK 1.2 Schmidt plates spanning 20 years was scanned 
with the COSMOS and SuperCOSMOS measuring ma- 
chines to provide a catologue of some 200,000 objects in 
the central 19 square degrees of the plate. These were cal- 
ibrated with CCD frames to provide light curves in Bj 
and R, and colours in UBVRI. Photometric errors have 
been discussed in detail in earlier papers (Hawkins 1996 
and references therein) and are ~ 0.08 magnitudes for any 
individual machine measurement, and only weakly depen- 
dent on magnitude. There are approximately four plates 
in each year which reduces the error to ~ 0.04 magnitudes 
per epoch. This is small compared with the amplitudes of 
interest in this paper, and no attempt has been made to 
deconvolve it. 



2.2. Sample selection 

For the analysis in this paper three samples will be 
defined. The first (UVX) is based only on position on 
the sky and ultraviolet excess. The area on the sky 
containing the sample is defined by a number of AAT 



AUTOFIB fields and a 2dF field (Folkes et al. |1999| and 
references therein), covering a total area of 7.0 square 
de grees. Within this area all objects with JJ — B < — 0.2 
and Bj < 21.0 were observed. This was extended to 
Bj < 21.5 in the 2dF field. The ultra-violet excess (UVX) 
cut although necessary to give a relatively clean sample of 
quasar candidates, has the well known limitation of only 
being effective for redshifts z < 2.2. Beyond this redshift 
quasars become red in U — B as the Lyman forest enters 
the U band. For samples at higher redshift, variability 
has been found to provide a very useful criterion for 
quasar selection (Hawkins 1996 ). The second sample for 
consideration herei (VAR), was selected on this basis, 
with the requirement that the object should lie in the 
same area of sky as for the first sample with a magnitude 
limit Bj < 21, or anywhere in the measured area of the 
plate with a magnitude limit Bj < 19.5, and should have 
an amplitude dm > 0.35. The defining epoch for the 
magnitudes was the year 1977. This sample was selected 
without any reference to colour and so can be used to 
measure trends over a large range of redshift (z < 3.5). 
The third sample (AMP) was selected in a similar 
manner, but over the whole measured area of the plate 
(19 square degrees), and with an amplitude cut dm > 1.1. 
There is some evidence that these large amplitude objects 
form a distinct group, which is discussed below. 

The development of fibre fed spectrographs has meant 
that every object included within a given set of search 
criteria can be observed to give a high completeness 



level. In the case of fields observed with AUTOFIB, 
the existence of forbidden regions in the 40 arcminute 
field, and the very variable throughput for different fibres 
aligned at different positions on the spectrograph slit 
meant that up to 20% of spectra did not have sufficient 
signal to provide an unambiguous redshift. These objects 
where re-observed later with the faint object spectrograph 
EFOSC on the ESO 3.6m telescope at La Silla. The 2dF 
observations where of very uniform quality, and the 
redshift measures or other classification were almost 
100% successful from the fibre-feed spectra. 

Details of all three samples, containing a total of 
384 quasars, are given in the Appendix. B magnitudes 
are in the Bj system defined by the Illa-J emulsion 
and the GG395 filter, and refer to the year 1977. The 
amplitude 8m is the difference between maximum and 
minmum magnitude achieved over the 20 year run of 
data. The samples which each quasar belongs to are 
indicated by 1, 2 and 3 corresponding to VAR, UVX and 
AMP respectively, as defined above. Data for many of 
the quasars have already been published by Hawkins & 
Veron p95| ), but are given again here for completeness. 
Any small differences in the parameters are the result of 
further refinement of the calibration, and more extended 
monitoring of the light curves. 

UVX selected samples have been used many times 
in the past for quasar surveys, and the constraints are 
quite well known. Variability selection has been less often 
used, and some additional comments are appropriate. 
In particular, there is the important question of com- 
pleteness. Hawkins & Veron ( 1995| ) use a small sample 
of quasars in field 287 from Morris et al. (1991) selected 
by objective prism to test completeness. In 1993, 79% 
of the objective prism quasars would have also been 
detected according to the variability criterion 5m > 0.35; 
by 1997 this figure had risen to 93%, with two quasars 
remaining below the variability threshold. In fact both 
of these quasars are clearly variable, with amplitudes 
0.32 and 0.29 magnitudes. Strangely, they lie at either 
extreme of the redshift and luminosity range of the 
sample, with redshifts 3.23, 0.29 and luminosities -27.55, 
-22.33 respectively. Fig. |l|(a) shows the distribution of 
epoch at which quasars first satisfy the detection criterion 
given by Hawkins ( |l996[ ). This is generally speaking 
equivalent to attaining an amplitude of 0.35 magnitudes. 
The distribution peaks between 2 and 4 years, and most 
quasars have satisfied the criterion after 7 years. Fig. |l|(b) 
shows the cumulative distribution of detection epochs, 
illustrating the point that after 15 years nearly all quasars 
have varied sufficiently to put them in the variability 
selected sample. 

The distribution of amplitudes is best shown with the 
UVX sample, which was selected without any reference 
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Fig. 1. The top panel shows a histogram of the epochs at 
which quasars could first have been detected on the basis 
of the variability criterion. The bottom panel shows the 
same plot in cumulative form. 
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Fig. 2. The top panel shows a histogram of amplitudes 
for the complete UVX detected sample. The bottom panel 
shows a similar histogram for the variability detected sam- 
ple (VAR). Quasars are constrained to have an amplitude 
greater than 0.35 mags. 



to variability. Fig. Q(a) is a histogram of amplitudes 
over a 20 year baseline, and shows a median amplitude 
of around 0.7 magnitudes. As would be expected from 
Fig. [l|, nearly all of the sample lies above the threshold 
amplitude of the variability selected sample of 0.35 
mags. With this in mind it is worth investigating the 
possibility of using the variability selected sample to look 
for correlations at higher redshift. Fig. ||(b) shows a 
histogram of amplitudes for the VAR sample, which from 
the way it was selected has a cut-off at an amplitude of 
0.35 mags. More interestingly, the distribution as a whole 
peaks at an amplitude of about 0.7, similar to the UVX 
sample. This peak is well clear of the cut-off, and implies 
that only a small fraction of quasars are missed from the 
VAR sample. There may however be a problem detecting 
the most luminous quasars, (Mb < —27) for which there 
is evidence for relatively small variations (Cristiani et al. 
1996). Nonetheless, the sample can be used with caution 



to extend quasar correlations to high redshift. 



3. Amplitude correlations 

The top two panels in Fig. || show the distribution of 
amplitudes as a function of redshift and luminosity for the 
UVX sample. The most striking thing about these figures 
is the marked drop in amplitude towards low redshift and 
luminosity. This is illustrated more clearly in the top two 
panels of Fig. ||, where the data is binned in intervals 
of 0.5 in redshift and unit absolute magnitude, and the 
mean plotted with y/~N error bars. The fall towards low 
redshift and luminosity is highly significanti (a 3-cr effect), 
although it is not clear from these data alone whether it is 
primarily a luminosity or redshift effect. The correlation 
coeffecient for the top left panel of Fig. || with z < 0.5 
is 0.46, and for the top right panel with Mb > —22 is 
-0.67, confirming the reality of the correlation. There is 
also some evidence for a decrease in amplitude for the 
most luminous quasars (Mb < —25), and an even weaker 
decline for high redshift objects. Correlation coefficients 
for all the data in the top two panels are 0.12 and -0.10 
for left and right panels respectively, emphasising the 
weakness of the effect. This is another manifestation 
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of the well known degeneracy between redshift and 
luminosity. This point will be investigated further below 
by dividing the data into sub-samples. 

The bottom two panels in Figs |^ and ^ show similar 
plots for the VAR sample. The greater redshift range 
allowed by variability selection still shows little trend of 
amplitude with redshift, but the decrease in amplitude 
for luminous quasars is shown to continue to greater 
luminosities, although the effect is inevitably lessened by 
the absence of quasars with Sm < 0.35. 

It has been mentioned above that there is a degener- 
acy between redshift and luminosity. This results from 
the fact that in a magnitude limited sample high redshift 
quasars tend to be high luminosity objects, and vice 
versa. Thus a trend with one parameter will be mimicked 
by a trend with the other, and the true relation will be 
hard to disentangle. This degeneracy can in principle be 
broken byi binning the data in redshift and luminosity, 
and the result of doing this is shown in Figs. || and |[ 
The VAR sample was used as it covers a wider range of 
luminosity and redshift, making binning feasible. The left 
hand panel shows amplitude as a function of luminosity 
in two redshift bins, z < 1.5 and z > 1.5. Data for the two 
redshift ranges overlap nicely, and clearly show a decrease 
in the relation between amplitude and luminosity. The 
right hand panel shows amplitude as a function of redshift 
in two luminosity bins. In this case there is some evidence 
for an increase in the relation at low redshift, but it 
is essentially flat beyond z = 0.5 for both luminosity 
ranges. It is however worth noting that the less luminous 
quasars have larger amplitudes, as expected from the 
data in the left hand panel. It would thus appear that 
while for Mb > — 25 amplitude does indeed decrease with 
luminosity, it does not change with redshift for z > 0.5. 

Examination of Fig. Ej, especially the bottom panels, 
suggests that there may be population of low luminosity 
and/or low redshift quasars distinct from the parent 
population. This is particularly evident in the bottom 
left hand panel of Fig. S, which is uniformly populated 
between amplitudes of 1.1 and 1.8 up to a redshift of 2 
at which point there is a sharp cut-off with no ampli- 
tudes greater than 1.1 at higher redshift. To investigate 
this population with better statistics all variables with 
Sm > 1.1 in the measured area of field 287 were observed 
on the 3.6m at La Silla to confirm their identification as 
quasars and measure redshifts. This became the AMP 
sample. Fig. |?j shows amplitude as a function of both 
redshift and luminosity, and it will be seen that there 
is indeed a cut-off in redshift at z ~ 2 and Mb ~ —25. 
It is clear that this cut-off must in fact be related 
to luminosity. If it were a redshift cut-off there is no 
reason why such objects should not be seen with greater 
luminosity. On the other hand if it were a luminosity 



cut-off, then this combined with a magnitude limit will 
indeed produce an effective cut-off in redshift. 



4. Discussion 

Attempts to measure correlations of amplitude (or some 
related variability parameter) with redshift and luminos- 
ity have a long history, which is summarised by Hawkins 
( 1996| ). Among recent work, a useful place to start is 
with the paper by Hook et al. ( 1994 ). They analyse a 
sample of ~ 300 quasars in the SGP area from 12 UK 
1.2m Schmidt plates taken in 5 separate yearly epochs 
spanning 16 years. They find a convincing anti-correlation 
between their variability parameter a v (a measure of 
variation about the mean) and luminosity, but a much 
weaker anti-correlation with redshift. They attribute 
this to the degeneracy between redshift and luminosity 
in their sample. The same sample is re-analysed by 



Cid Fcrnandes et al. (1996) who use a variability index 



related to variance, and also one related to the structure 
function. They concur with Hook et al. ( 1994| ) that there 
is an anti-correlation between their variability indices and 
luminosity, but claim a positive correlation with redshift. 
This effect is not apparent to the eye, but is interpreted 
as a variability-wavelength dependence rather than an 
intrinsic variability-redshift dependence. The net result 
in an un-binned sample cancels out with the luminosity 
variability relationship. Cid Fernandes et al.'s claim for 
a positive correlation between amplitude and redshift 
appears to be motivated at least in part by expectations 
arising from a paper by Di Clemente et al. ( 1996 ). This 
interesting paper examines the relation between their 
variability parameter S\ (an amplitude based on the 
structure function) and wavelength. Their sample is com- 
posed of PG quasars (Schmidt and Green 1983), which 
are mostly low redshift and relatively low luminosity 
objects. With the help of archival IUE observations they 
find that S\ decreases with wavelength. This effect can 
clearly be seen in the light curve from the intensive 
monitoring programme of the Seyfert galaxy NGC 5548 



(Clavel et al. 1991), which has a larger amplitude at 
shorter wavelength. Fig. || shows the relation between 
amplitude and wavelength taken from the light curves, 
with a best-fit quadratic curve. 

The relation between rest wavelength and amplitude 
is essentially equivalent to the relation between redshift 
and amplitude, where one is seeing progressively shorter 
wavelengths at higher redshift. This is illustrated in the 
top panel of Fig. || which shows the UVX sample with 
two curves superimposed. The solid line is converted 
from Fig. || and does not appear to follow the trend 
of the data, but the large scatter makes it hard to 
construct a convincing test. Nonetheless it suggests that 
any relation which holds for Seyfert galaxies might have 
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Fig. 3. The top two panels show plots of amplitude versus redshft and absolute magnitude for quasars selected solely on 
the basis of ultra-violet excess. The bottom two panels are similar plots for quasars selected on the basis of variability. 



to be modified fo r qua sars. The dotted line is from Di 
Clemente et al. ( |l996| ), and shows a very small effect, 
which does nonetheless follow the flat distribution of the 
data. 



The decrease of amplitude for the smallest redshift 
and luminosity seen in the top two panels of Fig. |] 
has a number of possible explanations, ft could be the 
effect of the underlying galaxy dominating any change in 
nuclear brightness for low luminosity objects, it could be 
a consequence of the small optical depth to microlensing 
at low redshift, or it could be a consequence of the 
wavelength dependence of variability (Cristiani et al. ?). 
The present dataset is not adequate to settle the question, 
which is best done by looking at luminous quasars at very 
low redshift and Seyfert galaxies at high redshift. 



All the plots in Figs || and || show a trend of 
decreasing amplitude towards higher redshift or more 
luminous objects. Although the trend as a function of 
luminosity is more marked, the old problem of degeneracy 
makes it hard to say for certain that it is a luminosity 
effect which is being observed. However, if we look at Fig. 
H where the data are binned in luminosity and redshift we 
see that while there is no significant trend of amplitude 
with redshift in either luminosity bin, there is a marked 
inverse correlation between amplitude and luminosity. 

The relation between amplitude and luminosity is in 
agreement with that found in earlier work (Hook et al. 



1994, Hawkins 1996 , Cristiani et al. 1996) and may well 



turn out to be a useful way of distinguishing between 
various schemes for quasar variability. The evidence for 
a constant amplitude with redshift is more debatable. It 
would appear to be consistent with the early claim of 
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Fig. 4. Plots based on the same data as for Fig. 3, but binned in redshift intervals of 0.5 and unit absolute magnitude 
intervals. The error bars show the uncertainty in the position of the mean. 



Hook et al. ( 1994 ) for a weak anti-correlation with redshift 
which they ascribed to degeneracy with luminosity. It is 
also in agreement with the results of Cristiani et al. in 
the observer's frame. When they correct their structure 
function to the quasar rest frame they inevitably imprint 
a positive correlation between their variability parameter 
and redshift. 

To investigate this dependence of amplitude on 
redshift for the present sample, Fig. [To] shows the epoch 
at which quasars achieve their maximum amplitude as 
a function of redshift. Apart from very low redshifts 
(z < 0.3) this relation is flat, implying that at least over 
the 21 years of the present dataset, time dilation effects 
will not bias the measurement of amplitude. Since the 
conclusions of Cid Fernandes et al. ( 1996| ) are largely 
based on a sample for which a time dilation correction 
has been applied, it is not feasible to make a direct 



comparison with the present work. However, it appears 
that the main difference between their results and those of 
Hook et al. is in the definition of a variability parameter 
and the method of analysis (both papers are based on the 
SGP sample). 



There are perhaps three currently discussed schemes 
for quasar variability. The least well constrained is the 
accretion disk model, where instabilities are propagated 
across the disk leading to variation in light. The de- 
tails of this approach have proved hard to work out, 
especially in the context of the constraints imposed by 
existing observations, but it does not seem to lead to an 
inverse correlation between amplitude and luminosity. 
An interesting recent attempt to model variation on the 
basis of accretion disk instabilities by Kawaguchi et al. 
( 1998 ) may provide a means for producing the observed 
variations. It does however appear to predict variations 
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Fig. 5. Plots of amplitude versus redshift and absolute magnitude for the variability detected sample (VAR). The 
top two panels show amplitude versus absolute magnitude for two redshift ranges, and the bottom two panels show 
amplitude versus redshift for two absolute magnitude ranges. 
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Fig. 8. Relation between amplitude and wavelength for 
the Seyfert galaxy NGC 5548 from HIE data published by 
Clavel et al. 1991. 



which are either too asymmetric or of too small an 
amplitude to be consistent with the current observations. 
The timescales which they predict are also rather short, 
around 200 days for reasonable input parameters, and 
much shorter than the observed timescale of a few years. 



An alternative approach, developed by Terlevich and 
his collaborators, accounts for the variation by postulat- 
ing that the quasar is powered by a series of supernova 
explosions. Qualitatively, this model can account for the 
observed relation between luminosity and amplitude, 
and works quite well for Seyfert galaxies (Aretxaga & 
Terlevich 1994). However, for quasars (Aretxaga et al. 



1997) large numbers of supernovae are required to achieve 
the luminosity, which results in smaller variations. For 
example, even a relatively modest quasar with absolute 
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Fig. 6. Plots based on the same data as for Fig. 5, but binned in unit absolute magnitude intervals and redshift 
intervals of 0.5. The left hand panel shows the relation between amplitude and absolute magnitude for quasars with 
z < 1.5 (small dots) and z > 1.5 (large dots). The right hand panel shows the relation between amplitude and redshift 
for quasars with Mg > —25 (small dots) and Mb < — 25 (large dots). The error bars show the uncertainty in the 
position of the mean. 
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Fig. 7. Plots of amplitude versus redshift and absolute magnitude for the most variable quasars in the sample. 



magnitude Mb ~ —26 would need some 300 type II 
supernovae per year to power it, which given typical 
decay times would lead to very little variation at all. 
This clearly conflicts with observations in this paper 
which show that most quasars vary by around 0.5 to 1 
magnitude on a timescale of a few years. 

The third way of explaining quasar variability is to 
invoke microlensing. This approach has been explored in 



several recent papers (Hawkins 1993, 1996, Hawkins & 



Taylor 1997), and seems to account well for a number of 



statistical properties of quasar light curves. It can also 
explain the inverse correlationi between luminosity and 



amplitude in a natural way. It is well known that when 
a point source is microlensed by a population of compact 
bodies of significant optical depth, the lenses combine non- 
linearly to form a caustic pattern which produces sharp 
spikes in the resulting light curves (Schneider & Weiss 
1987). As the source becomes comparable in size with, 
or larger than, the Einstein radius of the lenses the am- 
plitude decreases (Refsdal & Stabell 1991). Thus if one 



assumes a uniform temperature for quasar disks, and that 
the luminosity is determined by the disk area, the larger 
more luminous disks will be amplified less. Using a re- 
lation between source size (in terms of Einstein radius) 
and amplitude given by Refsdal & Stabell (1991) (eqn 1), 
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Fig. 10. Plot of the epoch at which the largest magnitude 
difference is achieved in quasar light curves at different 
redshifts. Poisson error bars are shown based on the num- 
ber of objects in each bin. 

does imply large amplitude variations for low luminosity 
objects, where each event is of comparable brightness to 
the nucleus itself. Perhaps the most natural explanation 
comes from microlcnsing, where the nucleus of low 
luminosity quasars would plausibly become very small 
compared with the Einstein radii of the lenses, resulting 
in large amplifications from caustic crossing events. 



Fig. 9. Plots of amplitude versus redshift and luminosity 
for UVX selected quasars. In the top panel the solid line 
represents the relation from Fig. 8, with wavelength con- 
verted to redshift; the dotted line is the model from Di 
Clemente et al. (1996). 



one can thus derive a relation between amplitude 5m and 
absolute magnitude M of the form: 

5m = 10°- 2 ( M+C ) 

where c is a constant. The bottom panel of Fig. ^| shows a 
plot of amplitude versus absolute magnitude for the UVX 
sample with this relation superimposed. The constant 
c was adjusted to allow the curve to track the upper 
envelope of the points, which has the effect of defining 
the quasar disk size. This ranges from 1.1 Einstein radii 
for M B = -23 to 11 for M B = -28. The points scatter 
downwards from the upper envelope because the quasars 
do not necessarily attain their maximum possible ampli- 
tude. Although one cannot say that the curve provides a 
fit to the data, the trend is certainly well represented. 

The possibility of a distinct population of large 
amplitude, low luminosity quasars suggested by Fig. ^| 
may also be used to test the models of quasar variability. 
Again, an accretion disk provides no obvious mechanism 
for such an effect. The Christmas Tree model certainly 



5. Conclusions 

In this paper we have examined the dependance of 
amplitude on redshift and luminosity for large samples of 
quasars selected on the basis of ultra-violet excess and 
variability. The quasars span a redshift range < z < 3.5 
and a luminosity range —20 > Mb > —28. There is 
evidence for a correlation of amplitude with luminosity 
and/or redshift for the sample as a whole, but when 
it is binned in redshift the correlation with luminosity 
becomes significant. This result could be strengthened by 
the possible existence of a population of large amplitude 
low luminosity objects in the sample. No convincing 
evidence is found for a correlation between amplitude 
and redshift, either for the sample as a whole or when it 
is binned in luminosity. 

Various models of quasar variability are examined 
with respect to the observed correlations. It is concluded 
that any straightforward interpretation of an accretion 
disk model is incompatible with the data. The Christmas 
Tree model has some merits, especially in the regime 
of low luminosity quasars and Seyfert galaxies, but for 
luminous quasars the rate of supernovae required is too 
large to be compatible with the observed variability. The 
microlensing model can be used to explain all the data, 
although it does require that the Einstein radius of the 
microlensing bodies is comparable in size to the quasar 



nucleus. 
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